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We present a computation of Cu K-edge resonant inelastic x-ray scattering (RIXS) spectra for 
electron-doped cuprates which includes coupling to bosonic fluctuations. Comparison with exper- 
iment over a wide range of energy and momentum transfers allows us to identify the signatures 
of three key normal-state energy scales: the pseudogap, charge transfer gap, and Mott gap. The 
calculations involve a three band Hubbard Hamiltonian based on Cu d^2_y2 and O Px, Py orbitals, 
with a self-energy correction which arises due to spin and charge fluctuations. Our theory repro- 
duces characteristic features e.g., gap collapse, large spectral weight broadening, and spectral weight 
transfer as a function of doping, as seen in experiments. 
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Cuprates are widely believed to be charge-transfer 
insulators Ij, with a Mott gap between Cu-d^2_y2 or- 
bitals much larger than the charge transfer gap between 
Cu-d and 0-p orbitals. The upper (UHB) and lower Hub- 
bard band (LHB) of the Cu orbitals are intimately related 
to the antibonding and bonding bands of the three band 
model, and it is important to understand how the strong 
correlations of Mott physics modify these bands from the 
conventional LDA-based picture. Meanwhile, a third en- 
ergy scale, the pseudogap scale, has been found exper- 
imentally, and its origin and relation to the other two 
scales continues to be a matter of intense debate. Here 
we model electron-doped Nd2-xCexCu04±5 (NCCO), for 
which the pseudogap is well described as a competing an- 
tiferromagnetic (AF) order. 

Experimental access to the LHB and/or the bonding 
band has proven difficult and the corresponding optical 
interband transitions have not been observed. Moreover, 
while the antibonding (ij.2_j^2 band lies at the top of the 
d-bands, the bonding dx2-y2 band is found in LDA to lie 
at the bottom of a veritable 'spaghetti' of c?-bands and 
their associated oxygen orbitals, nearly 6 eV below the 
Fermi level as seen in Fig. [T] Therefore, it is difficult 
to extract this band from ARPES data. On the other 
hand, RIXS is a local probe directly rearranging the Cu 
and O orbitals, and as such can provide selective access to 
the bonding bands. Indeed, RIXS experiments report a 
strong feature in most cuprates in the 6-8 eV range which 
has been associated with this band t2j-i7j. In this article 
we show that by incorporating strong renormalization of 
the near Fermi energy bands by magnon fluctuations [5] 
12], the high-energy RIXS features are indeed consistent 
with transitions from the LHB to the UHB. This resolves 
a puzzling discrepancy in earlier calculations [lOj which 
were unable to fit both the low and high energy parts of 
the spectrum. We also capture another important feature 
of the spectrum, the realistic broadening, which arises 



due to the strong coupling to bosonic quasiparticles. 

Remarkably, we find that all three energy scales are 
strongly influenced by the Hubbard U. The three energy 
scales are the following: 1) the Mott gap scale which is 
the result of transitions from LHB to UHB, 2) the charge 
transfer gap scale which persists as a residual feature into 
the overdoped regime and 3) the pseudogap or AF gap 
scale which collapses in a quantum critical point near op- 
timal doping. For convenience, we label the AF-split sub- 
bands of the antibonding bands as the lower (LMB) and 
upper ( UMB) magnetic bands. Cuprate magnetism nat- 
urally separates into two regimes: at high energies Mott 
physics produces localized spin singlets on each copper 
site, splitting the Cu dispersion by an energy ^ U into 
upper and lower Hubbard bands. In the presence of hy- 
bridization with oxygens, the LHB [UHB] becomes iden- 
tified with the bonding [antibonding] band of the three- 
band model. At lower energies, these singlets interact on 
different sites, leading to magnetic gaps in both UHB and 
LHB of magnitude ~ mj^U via more conventional Slater 
physics associated with AF order, where is the magne- 
tization on Cu. The Mott physics arises as an emergent 
phenomenon. When the AF gaps open at half filling, 
hybridization between Cu and O is mostly lost. For in- 
stance, in the antibonding band electrons in the UMB 
have mainly Cu character, while the opposite happens 
in the bonding band [TD]. Consequently, the states near 
the top of the lower magnetic band are of nearly pure 
oxygen character [10]. Thus, due to strong correlations, 
the 'charge-transfer' gap at half filling coincides with the 
AF gap. At finite doping, these two features separate in 
energy: the AF gap collapses rapidly pjj, while a resid- 
ual charge-transfer gap persists in optical spectra at high 
energies, due to strong magnetic fluctuations, closely re- 
lated to the high energy kink (HEK), or 'waterfall' effect 
seen in photoemission [H [12]. Here we show that this 
residual charge-transfer gap is also present in RIXS. 
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In K-edge RIXS the incident x-ray excites a Cu Is — >■ 
4p transition with an intermediate state shakeup involv- 
ing mainly Cu d^'^-y-^ and O p states. Within the RPA 
framework, the RIXS cross section for this process is 

[iniiiiiii] 

W(q,t^,w,) = {2nfN\w{u:,uj,)\'' 

Im[r+-(q,c.)] |aj2cos(2q.R4l) 

Here oji is the initial photon energy (taken to be -5 eV) 
and w, q are the energy and the momentum, respectively, 
which are transferred in the RIXS process. it;(a;,a;i) con- 
tains all the matrix-element information of the initial 
and final state transition probabilities [10^, N is the to- 
tal number of Cu atoms and is the position of the 
^th orbital present in the intermediate state. The near- 
est neighbor (NN) O excitations and second NN Cu ex- 
citations are included via ai and a2, respectively. We 
assume small values of ai = 0.1, = 0.05 in this 
study, whereas ao is equal to 1. The transferred mo- 
mentum and energy are then shared by the electron-hole 
pair created in the intermediate state from Cu dx2-y2 
and O p bands, Y^~{q,uj). In the Keldysh formal- 
ism it takes the form of a charge correlation function 
or the joint density of states (JDOS) (in real time) as 

- = (Pq,/xV'(t')P-q,M'T(0) ^Ith Pq,^<,(t) 

representing the charge operator. It is straightforward 
to show that Y can be calculated as the convolution be- 
tween the spectral weights (A) of the filled and empty 
states [131 [TS] as 

'^t!^a\f,ai^^^) ^ ^ j j duj2Af,i,,{k + q,UJi) 

k 

X A^,>^{k,UJ2) T^^T (2 

UJ + 10 + U)2 — ijJl 

where /(w) is Fermi function and cr is the spin index. The 
prime in the /c— summation means that the summation is 
restricted to the AF zone. 
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FIG. 1: (Color online) (a) Spectral weight of electronic states 
in NCCO for x — 0.14. A constant imaginary part of magni- 
tude 0.5 eV is added to broaden the spectra, (b) Calculated 
self energy for anti-bonding bands. The red dashed curve and 
the blue sohd curve are respectively real and imaginary parts 
of the self-energy, (c) Spectral weight as in (a), but modified 
by the self energy of (b). 



RIXS spectra are calculated using Eq. [2] in which the 
spectral weight A is computed using a three-band Hub- 
bard model with the Hamiltonian: 

3 i 

+ Yl tcuoid]p: + (cc)) + J2 tooipjp', + ic4^) 

<i,j> <i,i'> 

where A^q is the (bare) difference between the onsite 
energy levels of Cu dx2_y2 and O p — cr, tcuO the copper- 
d oxygen-p hopping parameter, too the oxygen-oxygen 
hopping parameter and U (Up) the Hubbard interaction 
parameter on Cu (O). Uj = d^^dj and rii = p\pi are 
the number operators for Cu-d and 0-p electrons, re- 
spectively. The equations were solved at Hartree-Fock 
(HF) level to obtain a self-consistent mean-field solution. 
Hartree corrections lead to a renormalized Cu-0 split- 
ting parameter A = A^o + Uud/I — UpUpjl^ where 
(rip) is the average electron density on Cu(0) [16]. The 
AF order splits the three bands into six bands as seen 
in Fig. [ija) . Since self-energy corrections are explicitly 
included, we use bare LDA-like dispersions [HI in the 
three-band model rather than the dressed, experimen- 
tal dispersions TO]. Thus hopping parameters are taken 
from LDA while interaction parameters U and A^g are 
adjusted to optimize agreement between the antibonding 
band splitting and earlier one band results [TUl fTMTO] . 
When this is done, we find that A^o is small and nega- 
tive while I] has a very weak doping dependence |20) . 

The renormalization of the antibonding band due to 
bosonic fluctuations is calculated via a self-energy based 
on the QP-GW formahsm [H], 



i;(k,iiLj„) = 



E 



dj 

2^ 



rG(k-q,zw„ + w')w^(q,w')- (4) 



Here F is the vertex correction and W ~ U^X is the in- 
teraction term which includes both spin and charge fluc- 
tuations. AF order in included along the lines of Ref. [T^] 
where the effective AF gap is kept the same as in the 
one band model. Finally the self-energy (E) is incorpo- 
rated into the three band dispersion via Dyson's equation 
G^^ ~ Gq^ — E and A is computed from the dressed 
G. Our calculation includes only fluctuations associated 
with the band closest to the Fermi level, which produces 
negligible broadening for w > 4 eV. Therefore for higher 
energy bands we include a constant broadening, S — 0.5 
eV, consistent with the ARPES data [22j . 

Figure [l] shows how self energy effects modify the 
dispersion of the various bands of the three-band AF 
model for x — 0.14, comparing bare (a) and dressed (c) 
bands. The imaginary part of the self-energy, plotted in 
Fig. [ijb), attains a maximum around 1.7 to 2 eV, which 
leads to a strong broadening of the spectral weight in this 



3 



energy range, both below and above the Fermi level (de- 
noted by the black line), producing a characteristic kink 
or 'waterfall' effect in the dispersion. We will see in con- 
nection with Fig. [2] below that this 'waterfall' effect leads 
to a significant broadening in the RIXS spectrum since 
the spectrum of Eq. [2] involves a convolution of the filled 
and empty states. Fig. [ijc) also shows that the self- 
energy softens the low energy bands nearest the Fermi 
level. This renormalization should also show up in the 
lowest branch of the RIXS spectrum, but this is restricted 
to very low energies and does not appear prominently in 
Fig. [2] 

Figure [2] shows the calculated RIXS spectra of NCCO 
for j: = and x = 0.14, reflecting the modulation of 
the spectral intensities of Fig. 1 via matrix element ef- 
fects, which are well known to be important in various 
highly resolved spectroscopies. H^^Hj Frames (a) and (d) 
include AF order but without self energy corrections, 
whereas the calculations in frames (b) and (e) include the 
self energy. The high intensities at energies around 5.6 
eV involve the transition from the lower Hubbard band 
to the unoccupied states of the antibonding band, reflect- 
ing the Mott gap feature. This '6 eV feature is present 
for all dopings. At half filling, in frames (a) and (b), the 
high intensities around 2 eV occur due to the transition 
within the antibonding Cu-0 band across the AF gap. 
This gap collapses with doping and as a result we find a 
smaller AF gap at 14% electron doping in frames (d) and 
(e), close to the QCP, consistent with earlier results [T4] . 
A key result is that the self energy produces a realistic 
broadening comparable to that observed experimentally. 
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FIG. 2: (Color online) RIXS spectra from NCCO for x ^ 0: 
(a) theory without and (b) with self energy corrections, and 
(c) experiment PJJ. (d)-(f) Similar figures for x — 0.14. 

In the RIXS calculations, the 6 eV feature is the most 
intense in the spectrum, consistent with most early ex- 
periments on a variety of cuprates [U , but more 
recent experimentsl31j. including those of Figs. 2(c,f)[52], 



employ a range of uji where the 6 eV feature is suppressed 
and the lower energy features can be more easily probed. 
Except for this feature, most features in the calculated 
RIXS intensities follow the experimental trends. In the 
undoped cuprate in Fig. [2] (b) , we observe a broad peak 
at F around 2.5 eV, with the intensity decreasing around 
the zone corner (tt, tt) while it remains strong around 
(tTjO). a similar level of agreement is found in the case 
of 14% electron doping in Fig. [2] in panels (d)-(f) [55] . 
The black dots in panels (b) and (e) represent the peaks 
of the experimental spectra, reproducing the blue, black, 
and purple dots of panels (c) and (f). The agreement is 
remarkable for both dopings, a; = in frames (b) and (c) 
and X = 0.14 in frames (e) and (f) [34]. Results for 
a;=0.09 are similar, and are omitted for brevity. 

We comment here on the three energy scales. While 
the dispersions which follow from Eq. 3 are rather compli- 
cated, we find numerically that the Mott gap is approxi- 
mately equal to U and the AF gap to Uirid, as illustrated 
by arrows in Fig. 2(a). Also, the charge transfer energy 
is the difference between the average oxygen energy and 
the upper Cu band [T], which we find to be ~ U/2. Thus 
all three energy scales are controlled by U. In our calcu- 
lation the 6 eV feature represents transitions across the 
true Mott gap, and the good agreement with experiment 
indicates that RIXS can be used to probe this important 
scale and how it is modified by hybridization with oxy- 
gens - is the bonding band split as our calculations sug- 
gest? This feature will be discussed further below when 
we describe fits to individual g-cuts of the RIXS spectra. 
In optical spectra at half filling the ~2 eV charge trans- 
fer gap is indistinguishable from the AF gap[5]. At finite 
doping these two features separate, with the AF gap re- 
flected as a midinfrared peak which collapses rapidly with 
doping, while a residual charge transfer gap persists as 
a weak feature near 2 eV in the strongly doped regime. 
A similar evolution is found in RIXS. The RIXS lead- 
ing edge follows the doping dependence of the AF gap 
[TUl [T^ , while in Fig. 3 we show that a residual charge 
transfer gap feature can be seen in the RIXS spectra near 
the F point. Our three band calculation successfully re- 
produces the experimental finding that the magnetiza- 
tion scales with the AF gap [321 ISS]- Our calculated 
three-band RIXS spectrum modified by self-energy beau- 
tifully displays the broad feature around 2 eV at F in 
panel (b), in good agreement with experimental results 
in panel (a). Panels (c) and (d) display RIXS intensi- 
ties obtained from theory (blue) as well as experiment 
(red dots) as a function of energy at F, compared with 
the optical spectra [37] (green dashes) for x = 0.10 and 
a; = 0, respectively. The peak of experimental intensity is 
shifted towards slightly higher energy than the theoret- 
ical intensity, but the broadening is comparable in both 
cases. 

For more quantitative estimates of the broadening. Fig. 
Hcompares theoretical (blue solid) and experimental (red 
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FIG. 3: (Color online) RIXS spectra from NCCO for x = 
0.075, (a) experiment [29] and (b) theory. Intensity cuts along 
r, (c) for X = 0.075 and (d) for a: = with RIXS theory 
(blue solid line), RIXS experiment (red dot-dashed line) and 
optical experiment (green dashed line). In frames (a) and (b), 
solid arrows indicate intensity peaks along F, while dashed 
arrows indicate cuts taken in frame (c). In frames (c) and (d), 
the theoretical curves have been convoluted with a Gaussian 
broadening of 200 meV to mimic experimental resolution. 
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FIG. 4: (Color online) Comparison of theoretical (blue solid 
line) and experimental (red dashed line) spectra at half- 
filling with q fixed at (a): (7r,0), (b): F, (c): (7r/2,0), (d): 
(7r/2,7r/2). Panel (b) displays the high-energy RIXS peak 
from Ref. [S] The theoretical curves have been convoluted 
with a Gaussian broadening of 200 meV to mimic experimen- 
tal resolution. 



dashed line) RIXS intensities as a function of oj for several 
constant g-cuts. There is an overall good agreement in 
peak positions as well as lineshapes and broadening for 
all momenta. In particular, panel (b) shows the high 
energy RIXS feature from Ref. [2j. The good agreement 
with theory strongly suggests the identification of this 
feature with the LHB in NCCO. A similar peak is found 
in all cuprates, as would be expected for Mott physics. 

In conclusion, we find that RIXS is a suitable probe 
across all energy scales, including AF gap, charge- 
transfer, and Mott physics. We provide a three-band 
model that is capable of explaining the experimental 
RIXS spectra over the entire energy and doping range. 
We find a good correspondence between the RIXS spec- 
tra at F and the optical spectrum, but RIXS has the 
additional advantage of full momentum-space resolu- 
tion. While we have concentrated on the electron doped 
cuprates, our model should apply equally well to the hole 
doped case. 
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